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Abstract 


To account for surface reflection is a major and most complex issue of any satellite aerosol 
retrieval algorithm. The surface contribution to the satellite signal is usually calculated assuming 
Lambertian reflection. Then the problem is reduced to the calculation of total atmospheric 
transmittance and spherical albedo. Usually, these parameters are stored in so-called look-up tables 
(LUTs), which can be quite large. The main aim of this paper is the parameterization of these look-up 
tables in terms of simple approximate equations, which can be subsequently used in aerosol remote 
sensing techniques over reflecting land surfaces. 
© 2004 Elsevier B.V. All rights reserved. 
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1. Introduction 


The accuracy of aerosol remote sensing techniques heavily depends on the information 
on the ground reflectance. It is usually assumed that the ground can be approximated by a 
Lambertian surface, which simplifies the calculations to large extent. In particular, the 


* Corresponding author. Tel.: + 49 421 218 4475; fax: +49 421 218 4555. 
E-mail address: alexk@iup.physik.uni-bremen.de (A.A. Kokhanovsky). 


0169-8095/$ - see front matter © 2004 Elsevier B.V. All rights reserved. 
doi: 10.1016/j.atmosres.2004.07.004 


38 A.A. Kokhanovsky et al. / Atmospheric Research 73 (2005) 37—43 


account for the Lambertian ground reflection in aerosol remote sensing problems can be 
done using standard equations (Sobolev, 1972): 

Atib | T=T+ At;|R (1) 
1-rA 1-rA 
valid for homogeneous light scattering media. Here R is the aerosol reflection function, T 
is the aerosol transmission function, R,=R(A=0), 7;,=7(A=0), A is the ground albedo, r is 
the aerosol spherical albedo, ¢, is the total (diffuse+direct) atmospheric transmittance from 
the sun to the surface and f, is the total transmittance from the surface to the satellite. R is 
the diffuse reflectance in the direction of the satellite which approaches R =1—t for non- 
absorbing aerosol. The exact definition of all parameters and functions in Eq. (1) is given 
by Sobolev (1972). 

Usually R, and T; are tabulated (after subtraction of the Rayleigh contribution). Then the 
aerosol optical thickness t can be retrieved from the look-up tables for a given observation 
geometry, and prescribed aerosol phase function and single scattering albedo œo. It is often 
assumed that the single scattering albedo is close to 1. For a non-zero surface albedo, one has 
to account for the contribution of the surface reflection, which is given by the second terms in 
Eq. (1). The value of A is not known; hence, an a-priori assumption of the spectral surface 
albedo A(A) is required. There are also possibilities to tune this assumed spectrum as 
proposed by Von Hoyningen-Huene et al. (2003). In particular, tı, t2, r, and R are 
tabulated for various observation/sun geometries, aerosol single scattering albedos, optical 
thicknesses and phase functions. This leads to rather large, multi-dimensional databases. 
The aim of this paper is the reduction of these look-up tables (LUTs) by developing a 
parameterization based on exact radiative transfer calculations. Such a parameterization 
allows a fast online calculation of t), t2, r, and R for exactly the required conditions— 
avoiding the need to pre-calculate extensive databases and to interpolate between the 
stored grid points. Only the case of nonabsorbing media is considered. The uncertainty of 
the proposed parameterization is well below the uncertainty of the spectral surface albedo 
A(2). Hence the error is negligible as far as practical aspects of passive aerosol remote 
sensing from space are of concern. The equations can also be used for back-of-envelope 
calculations of the surface influence on atmospheric reflectance and transmittance. 


R= R, + 


2. The parameterization of total transmittance and spherical albedo 


In a first step, diffuse transmittance and spherical albedo were calculated with the 
libRadtran radiative transfer package developed by Kylling and Mayer (1993-2004). 
LibRadtran includes a detailed atmospheric model, with molecules, aerosols, as well as 
water and ice clouds, which can be defined and modified in a user-friendly and flexible 
way. The package provides different numerical solvers of the radiative transfer equation, 
including exact discrete-ordinate schemes, a fast two-stream code, a polarization- 
dependent solver, and even a fully three-dimensional Monte Carlo code for horizontally 
inhomogeneous problems. LibRadtran has been extensively validated by comparison with 
other radiative transfer schemes (e.g. Van Weele et al., 2000) and by comparison with 
observations (e.g. Mayer et al., 1997). For the particular application presented here, the 
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discrete ordinate solver by Stamnes et al. (1988), version 2.0, was chosen which basically 
allows exact calculation of radiances and irradiances. The results obtained were compared 
with those derived from the SCIATRAN (Rozanov et al., 2002) radiative transfer code. 
The differences were found to be smaller than 1%, which confirms the high accuracy of 
both radiative transfer solvers. 

The commonly used Henyey—Greenstein scattering phase function was chosen for the 
aerosol which is an analytical function of the asymmetry parameter g: 


1—2? 
0) = : 2 
P(9) (1 + g2 — 2gcos0)!* (2) 


Here 0 is the scattering angle. Eq. (2) is only a rough approximation of a real aerosol 
scattering phase function but this is only of minor importance for our application because 
diffuse transmittance and reflectance are mostly governed by the value of the asymmetry 
parameter g and not by the precise angular structure of the phase function. 

We propose the following parameterizations based on the exact radiative transfer 
calculations performed for w=: 


r = t(aexp( — t/a) + bexp( — t/f) +c), (3) 


t = exp( — 1/u;) + texp( — uj — vt — wt’), (4) 


The second term in Eq. (4) gives the diffuse transmittance. Here u, is the cosine of the 
solar angle, u2 is the cosine of the satellite angle, and 


3 
uj = Shalt, (5) 
m=0 
vj = po + piexp( — p24), (6) 
wj = qo + qiexp( — q24;). (7) 


The constants a, b, c, &, P, Po, qo, P1» q1» P25 42, ho, h1, h2, h3 are parameterized using 
polynomial expansions with respect to the asymmetry parameter g of the Henyey— 
Greenstein phase function. The standard regression routine incorporated in the ORIGIN 
package has been used for the parameterization. For all constants a third order polynomial 
is used (except for ho where a fourth order polynomial is required), e.g. 


3 
a= > asg? (8) 
s=0 

for the coefficient a. Coefficients for all parameters are given in Table 1. Note that Eqs. (3) 
and (4) provide the correct behavior for r, t; in the limit t—0: In this case r—~0 and t;—i. 
The parameterization is valid for g<[0,0.9], u;<[0.2, 1.0], and t<[0,2]. This covers the 
usual variability of aerosol optical thickness and observation geometries. The asymmetry 
parameter of aerosol media is typically close to 0.7 (Kokhanovsky, 2004). Our technique 
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Table 1 

List of constants for Eqs. (3)-(8) (Ag4=—2.3) 

s 0 1 2 3 

ds 0.18016 —0.18229 0.15535 —0.14223 
bs 0.58331 —0.50662 —0.09012 0.0207 
as 0.16775 —0.06969 0.08093 —0.08903 
Bs 1.09188 0.08994 0.49647 —0.75218 
C3 0.21475 —0.1 0.13639 —0.21948 
hos —1.88227 0.53661 —1.8047 3.26348 
his 5.97763 —2.04621 —2.0173 1.44843 
hos —5.47825 2.42154 —3.37057 6.13805 
has 2.07593 —2.03761 6.25975 —7.35503 
Pos 0.4923 1.0471 —2.61112 1.53155 
Pis 4.01521 —0.25886 —2.85378 3.61515 
Pas 3.76447 3.29106 —12.37951 9.85 
dos 0.000076 —0.316 0.67744 —0.4093 
qis —1.31136 —0.8901 3.55 —3.0646 
42s 5.21931 7.2255 —23.43878 17.65629 


can be applied not only to aerosol media but also to thin clouds (g = 0.75—0.85) and to the 
case of isotropic scattering (g=0). Parameterizations for optically thick media are given in 
a separate publication (Kokhanovsky and Mayer, 2003). For most cases, errors of the 
parameterization are not larger than 2%, which is satisfactory for most aerosol remote 
sensing problems (especially taking into account other uncertainties involved). A detailed 
study of the accuracy of Eqs. (3) and (4) is given in the next section. 


3. The accuracy of the parameterization 


Let us consider the accuracy of the proposed parameterization for the spherical albedo r 
(see Fig. la). The spherical albedo strongly depends on the asymmetry parameter g. E.g., it 
can reach 0.15 for molecular scattering at t=0.2. It is around 0.35 for heavy aerosol loads 
with t=2.0 and g=0.7. The error of Eq. (3) for the spherical albedo r is rather small, see 
Fig. 1b. In particular, it is smaller than 2% for teE[0.01, 2.0] and ge[0.0, 0.9]. A 
parameterization of the spherical albedo r for larger values of t is given by Kokhanovsky 
and Mayer (2003). 

The total transmittance t is given in Fig. 2a as function of t at g=0.7 and several values 
of u. The error in the total transmittance ¢ is presented in Fig. 2b. We see that the errors are 
smaller than 2% for most cases. In particular, for g=0.7, it is smaller than 3% for u>0.3. 
Note that the vertical spread in Fig. 2b is due to different angles and that in Fig. 2c is due 
to various optical thicknesses from 0 to 2.0. We conclude that the maximum error is 
around 8% for g—[0.0, 0.9], t=[0,2], we=[0.2, 1.0]. The error is smaller than 4% for t<1.6, 
g<0.8, we[0.2, 1.0]. This is the most important range of parameters as far as aerosol 
remote sensing is concerned. 

Because of the relation R =1—f valid for nonabsorbing media, the errors of our 
parameterization for R can be found from those given in Fig. 2b for the value of t2. In 
particular, we have: dg = — Kô,, where K = b /R, ôr is the relative error in the value of 
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Fig. 1. (a) The spherical albedo as function of optical thickness for the Henyey—Greenstein phase function with 


various values of g=0.0—0.9. Crosses give the results of our parameterization, circles are obtained by exact 
radiative transfer calculations with @o=1. (b) The relative error of Eq. (3) in % for various values of g. 


R, and ô, is the relative error in the value of t,. We found that dz is generally smaller than 
5% (except, at u=l, where R is rather small and, therefore, K = 3.0|ôr| = 10%). 
The relative error ôr in the reflection function R (see Eq. (1)) is given as follows: 


Ör = My/ 5; + ò, + N°87, (9) 


where M=R,/R, R,=Att2/(1—rA), N=rA/(1—rA). The error decreases with A as one might 
expect. The relative error in the transmission function T can be determined with a similar 
expression. Then one needs to substitute ft. by R in Eq. (9). 


4. Conclusion 


Simple analytical equations for the spherical albedo, diffuse and total transmittance of 
aerosol media are given. A parameterization has been developed based on exact radiative 
transfer calculations for a homogeneous nonabsorbing layer characterized by the Henyey— 
Greenstein phase function. 

The parameterization can be used in aerosol remote sensing algorithms (Von 
Hoyningen-Huene et al., 2003). It can be applied also to related atmospheric research 
problems (e.g. for ocean-color algorithms; see e.g. Santer and Martiny, 2003). The 
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Fig. 2. (a) The total transmittance as function of the optical thickness for the Henyey—Greenstein scattering phase function with g=0.7, wo=1, and various values of u=0.2— 
1.0. Crosses give the results of the parameterization, circles are obtained from exact radiative transfer calculations for nonabsorbing media. (b) The relative error of Eq. (4) 
in % for various values of g=0.0—0.9 and w=0.2—1.0. The vertical spread of data is due to various angles u. (c) The relative error of Eq. (4) in % at g=0.7 as function of the 
cosine of the observation angle u. The vertical spread is due to the variation of the error with optical thickness. The largest errors occur for large solar and satellite 
observation angles. Such cases are rarely used in aerosol remote sensing problems. 
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uncertainty of the parameterization is typically below 2%, which is accurate enough for 
most of applications. 
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